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1 Introduction 



Quark masses are fundamental parameters of the Standard Model (SM) and thus it is 
desirable to determine their numerical values with the highest possible precision. In 
order to do so it is necessary to fix a renormalization scheme which defines the quark 
mass. A renormalization scheme which is of particular importance in those situations 
where the quark mass is small as compared to the typical energy scale is the MS scheme. 
Although it has no immediate physical interpretation, in general a good convergence in 
perturbation theory is observed. On the other hand, the most intuitive renormalization 
scheme is the on-shell scheme where the renormalized quark mass is defined as the pole 
of the propagator. It is well known that such a definition is plagued by numerically large 
long-distance effects which are connected to the renormalon corrections. As a consequence 
one observes a bad behaviour of the perturbative expansion if the pole mass is used as a 
parameter. Still, it plays a crucial role in all computations of threshold phenomena like 
energy levels or the decay of a quark- ant i- quark bound state. In order to make contact 
between high-energy phenomena on one side and threshold processes on the other side it 
is necessary to have a precise relation between the pole and the MS quark mass at hand. 
Important three-loop contributions are provided in this paper. As a by-product we also 
compute the same kind of corrections to the wave function renormalization constant. 

The one-loop corrections to the relation between the MS and on-shell mass have been 
computed almost 30 years ago in Ref . [1] . About ten years later, the two-loop corrections 
have been calculated in Ref. [2]. Shortly afterwards also the two-loop result for the on- 
shell wave function renormalization constant Z® s has been obtained [3]. Contributions 
with a single mass scale were expressed via master integrals exactly in d = 4 — 2e; the 
e-expansion of the only non-trivial master integral is known [4, 5] far enough to obtain 
0(e 4 ) terms in Z^f and Z° s . Contributions with two mass scales were obtained up to 
0(1) only. Later, a reduction algorithm for two-loop two-scale on-shell integrals has been 
constructed, and expressions via master integrals exact in d have been obtained [6]Q The 
three-loop relation between the MS and on-shell mass has been computed at the end of 
the nineties in a semi-numerical way in Refs. [9, 10] where the off-shell fermion propagator 
has been considered for small and large external momenta. The on-shell quantities have 
been obtained with the help of a conformal mapping and Pade approximation. Later 
the results have been confirmed in Ref. [11] by an analytical on-shell calculation. The 
three-loop result for Z® has been obtained in Ref. [12]. Both analytical calculations have 
recently been rederived in Ref. [13]. In this paper we complete the three-loop results for the 
mass counterterm and Z® s by computing the contributions where a second mass scale 
is present through a closed quark loop. We would like to mention that the approximation 
linear in the mass ratio has been derived in Ref. [14] using the corresponding corrections 
to the static potential [15] and their cancellation in the relation between the MS and IS 

1 The e-expansion of the two non-trivial master integrals was only obtained up to 0(1). One of these 
integrals was later expanded to 0(e 5 ) [7] (though only the 0(e) term is published). Both integrals were 
calculated [8] up to 0(e), as series in the mass ratio up to the sixth order. 
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quark mass. 

We concentrate on the situation where the second quark mass is smaller than the external 
one, although our formulae can also be applied to the reversed situation. Then, however, 
often it is advantageous to perform a decoupling of the heavy mass leading to an effective 
theory where the latter is integrated out. 

As far as the quark masses in the SM are concerned the corrections considered in this 
paper are of practical relevance for the bottom quark where the second mass scale is 
given by the charm quark. In this case we have for the mass ratio m c /mt> « 0.3. Thus 
the massless approximation does not provide a good result. In all other cases the effect 
from a lighter quark mass is negligible. Nevertheless, below generic results are presented. 

In principle at three-loop order there is also a contribution involving two additional masses 
which are present in two closed fermion loops. However, for all practical applications the 
lightest quark mass can safely be neglected. 

The remainder of the paper is organized as follows: In Section [2] we discuss the framework 
which is used for the calculation. Afterwards the results for the relation between the on- 
shell and MS quark mass and the wave function renormalization constant are discussed 
in Sections [3] and HI respectively. Since the analytical expressions are quite involved we 
present in both cases handy approximation formulae of our results. Finally, Section [5] 
contains a simple application and our conclusions. Appendix |A] lists all master integrals 
in graphical from. 



2 On-shell renormalization of quark mass and wave 
function 

The formulae relevant for the computation of the renormalization constants for the mass 
and wave function have been derived in Refs. [12,13]. For completeness we repeat the 
resulting expressions which read 

= 1 + M fl ), (1) 



(Z°Y = 1 + 2MJ— S 1 (g 2 ,M ? ) ?2=M2 +E 2 (M 9 2 ) M ? ), (2) 

where Z° and Z® are defined through 

m ?j0 = Z° S _M 9 , (3) 

'OS 



= y/Z™il>. (4) 

ip is the quark field renormalized in the on-shell scheme with mass m q , M q is the on-shell 
mass and bare quantities are denoted by a subscript 0. E denotes the quark self-energy 
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contributions which can be decomposed as 

T,{q,m q ) = m,£i(g 2 ,m 9 ) + (jtf - m,) S 2 (g 2 ,m g ) . (5) 
For completeness let us also introduce the MS renormalization constant via 

m Q,o = Z^ s m q {fi). (6) 

The quantities on the right-hand side of Eqs. ([1]) and (j2J) are obtained by considering the 
external momentum of the quarks to be q = Q(l + t) with Q 2 = M 2 . The application of 
the projector + M (? )/(4M ( 2 ) and an expansion to first order in t leads to 

Tr{^^S(g,M 9 )} = E^ 2 , M q ) + tE 2 (g 2 , M q ) 

= Ei(Mg,M g ) + (2M 9 2 A El(g 2 )M?) |^ +E2(M 2 )M ^ t 
+0(t 2 ) . (7) 

Thus, to obtain Z,^ s one only needs to calculate E^ for g 2 = M 2 . To calculate , one 
has to compute the first derivative of the self-energy diagrams. The mass renormalization 
is taken into account iteratively by calculating one- and two-loop diagrams with zero- 
momentum insertions. 

The results for the renormalization constants can be cast into the following form 
= l + ^( e —YsZ^+(^)U e —Y £ 5Z^ 
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with i e {m,2}. It is convenient to further decompose the two- and three-loop contribu- 
tion in terms of the different colour factors 

bzf ] (x) = C 2 F Z FF + C F C A Z FA + +C F T F mZf L + C F T F n h Z FH + C F T F n m Z FM {x) 
8Z® (x) = Cl Z FFF + C 2 F C A Z FFA + C F C\ Z FAA + C F T F m (C F Z FFL + C A Z FAL 
+T F n t Z FLL + T F n h Z FHL + T F n m Z FML (x)) 

+ C F T F n h (C F Z FFH + C A Z FAH + T F n h Z FHH + T F n m Z FMH {x)) 

+ C F T F n m (C F Z FFM (x) + C A Z FAM (x) + T F n m Z FMM (x)) , (9) 

where n\ and mark the closed quark loops with mass zero and m q , respectively. n m 
labels the closed quark loops involving a second mass scale which we denote as rrif. 
Although we have rih — 1 and n m = 1 in our applications, we keep a generic label 
which is useful when tracing the origin of the individual contributions. In Eq. ([9]) C F = 
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Figure 1: All three-loop heavy-quark self-energy diagrams containing n m . Thick solid 
lines denote massive quarks with mass m q and thin ones quarks with mass m,f. Wavy 
lines denote gluons and crosses mark counterterm insertions. The shaded blob denotes 
the sum of one-loop massless self-energy insertions (massless quarks, gluons, ghosts). 

— 1)/(2N C ) and Ca = N c are the eigenvalues of the quadratic Casimir operators 
of the fundamental and adjoint representation of SU(N C ), respectively. In the case of 
QCD we have N c = 3. Tp = 1/2 is the index of the fundamental representation and 
rif = rii + rih + n m is the number of quark flavours. a s (fj) is the strong coupling constant 
defined in the MS scheme with rif active flavours. The coefficients proportional to n m are 
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functions of the mass ratio which we define by 



Mr , s 

x = 4' (10) 



i.e. the ratio of the on-shell masses. The dependence of the coefficients in Eq. (Q on L 



log// /Mr is suppressed. The Feynman diagrams producing n m -dependent contributions 
to Z° and are shown in Fig. [IJ 

The one-loop result for expanded up to order e 2 can be found in Eq. (13) of Ref. [13] 
and the expression for 5Zm^\ nm =o including 0(e) terms is given in Eq. (14) of the same 
reference. Z^ M and Z% have been computed in analytic form in Ref. [2]. The main 
result of this paper are the functions Z™, Z FAM , Z* LM , Z* HM and Z™ M which are 
discussed in Section [3j 

In the case of the mass renormalization it is convenient to consider the ratio between the 
on-shell and MS renormalization constants 

z ms M„ 1 J 

which is finite. We furthermore adopt the notation introduced in Eqs. (jSJ) and (Q for 
Z° s also for z m . Let us for later reference provide already here the result for the MS 
renormalization constant which is given by [16] 

z™ = i+y>(^#Y, (12) 



i+£a(- 

i=l V - 



7T 



with 



1/15 1 \ 1 / 101 



1 / 65 7 1 2 \ 1 /2329 25 5 2 \ 

^ = ?l-l^ + ^-^J + ?lW-^ + 648^J 

1 / 1249 5 , 277 35 2 \ 

+ - 1 Gn f H 71* H n 2 f . (13) 

e V 192 18 s ; 648 ; 3888 V v ; 

As far as the wave function renormalization constant is concerned, we have at the one-loop 
level SZ^ = SZm and the two- loop contributions for n m = 0, including order e terms, 
can be found in Eq. (25) of Ref. [13]. The results for the functions Zg FM , Z FAM , Z FLM , 
Z FHM and Z F MM are discussed in Section H 



Starting from the three-loop level, the wave function renormalization constant depends 
on the gauge parameter, £, which we define through the gluon propagator as 

i I £ k^k,, 



where a and b are colour indices. 
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3 On-shell mass relation 



3.1 Two-loop result 

Before discussing in detail the three-loop result let us consider the two-loop quantity z™ . 
The analytical result can be found in Ref. [2] and is given by 

z™ = ^{48x 4 log 2 (x)+48x 2 log(x) + 72x 2 + 4L M (3L M +13) 
+ 8tt 2 (x 4 - 3x 3 - 3x + 1) + 71 

- 48(a; + l) 2 (x 2 - x + l) [log(s) log(x + 1) + Li 2 (-x)} 

- A8(x - l) 2 (x 2 + x + 1) [log(x) log(l -x) + Li 2 (x)] } , (15) 

where L M = \og(fi 2 /M 2 ) and Li 2 is the dilogarithm. In our approach all occuring inte- 
grals are reduced to four master integrals [6]. The Harmonic Polylogarithms (HPL) [17] 
which appear in a first step in the results of these integrals can be transformed into the 
(di)logarithms of Eq. (f!5l) resulting in complete agreement with Ref. [2] a 

In Fig. is shown as a function of x for /i 2 = M 2 . In addition to the exact result we 

also show the curves corresponding to the linear approximations for x — > and x — > 1. 
One observes a rapid convergence of both expansions which almost extends to x = 1 
(x = 0) once the x 10 ((1 — a;) 10 ) terms are included. Note that the x 5 terms are sufficient 
in order to provide an excellent approximation far below the per mill level for x = 0.3. 
We also want to mention that a 12% deviation is observed (for x = 0.3) if only the linear 
terms in x are included into the expansion. 

Let us finally provide an approximation formula [2] which agrees to better than 1% with 
the analytical formula of Eq. ffTS]) : 

13 1 

z™ = 1.562 - 2.394x + 0.9428x 2 - 0.2666x 3 + — L M + -L 2 t . (16) 

24 8 

3.2 Three-loop corrections 

At three-loop order we want to discuss in a first step the coefficients of the five ra m - 
dependent colour structures but afterwards also consider the physical quantity which is 
obtained after inserting the numerical values for the QCD colour factors. 

All Feynman diagrams are generated with QGRAF [18] and the various topologies are identi- 
fied with the help of q2e and exp [19,20]. In a second calculation the three-loop diagrams 
were generated starting from three generic topologies which already appear in the cal- 
culation of Ref. [21]. In a next step the reduction of the various functions to so-called 
master integrals (Mi's) has to be achieved. For this step we use the so-called Laporta 

2 A trivial overall factor Cf — 4/3 is actually missing in Eqs. (17) and (20) of this reference. 
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Figure 2: The two- loop correction z™ as a function of x at /i = M 9 . Next to the exact 
result (solid line) also the approximations for x — > (dash-dotted line) and for a: — ► 1 
(dashed line) are shown including the linear terms. 

method [22, 23] which reduces the three-loop integrals to 27 Mi's. We use the implemen- 
tation of Laporta's algorithm in the program Crusher [24] . It is written in C++ and uses 
GiNaC [25] for simple manipulations like taking derivatives of polynomial quantities. In 
the practical implementation of the Laporta algorithm one of the most time-consuming 
operations is the simplification of the coefficients appearing in front of the individual in- 
tegrals. This task is performed with the help of Fermat [26] where a special interface has 
been used (see Ref. [27]). The main features of the implementation are the automated 
generation of the integration- by-parts (IBP) identities [28], a complete symmetrization of 
the diagrams and the possibility to use multiprocessor environments. 

In Figs. [7HT01 of Appendix |X] a graphical representation of the master integrals can be 
found. We have chosen two independent ways to compute the e-expansion of the master 
integrals. The first one relies on the Mellin-Barnes technique (see, e.g., Ref. [29]) and 
provides us with numerical results. Here we have used the Mathematica package MB .m [30]. 
With the help of our second method, based on differential equations, Ref. [31], we were 
able to evaluate all but four master integrals in analytic form. More details can be found 
in Ref. [32]. 

The coefficient functions of all master integrals contain a 1/e pole, some even a 1/e 2 
term which means that the master integrals have to be expanded up to order e and e 2 , 
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respectively. All but four master integrals could be evaluated analytically in terms of HPLs 
which we evaluate numerically with the help of the Mathematica package HPL.m [33,34]. 
Two of the remaining four integrals, which are all needed up to order e, could be computed 
including the constant term [32] and for the residual two integrals analytical results are 
obtained for the pole parts. For the still remaining six coefficients integral representations 
are available which in the worst case are two-dimensional. 

Close to x = we observe large cancellations between the contributions originating from 
different master integrals. On the other hand, as we have seen above, the expansion for 
i < 1 converges very fast; at two- loop order the first five terms approximate the exact 
result to 0.02% for x = 0.3 which is relevant for the charm mass effects to the bottom 
quark mass. For this reason we decided to derive an expansion of our result including 
terms of order x 8 . The results expressed in terms of HPLs can simply be expanded using 
HPL.m [33,34]. For the remaining coefficients we use their Mellin-Barnes representation 
in order to express them in terms of multiple sums which in turn leads to the coefficients 
of x n . 

Our result for the five functions z^ FM , z^ AM , z^ LM , z^ HM and z^ MM are shown in Fig. [3] 
for < x < 1. The exact results are represented by thick lines and the expansion terms 
for x — ► up to the linear term as thin lines. The latter provide a good approximation to 
the exact results up to about x ~ 0.1 . . . 0.4, depending on the colour structure. We want 
to mention that the expansion terms including corrections of order x 8 provide a good ap- 
proximation almost up to x = 1. Note that for x = 0.3 the linear approximation deviates 
from the exact result (obtained by the proper sum of the individual colour structures) by 
8%, whereas the deviation including terms up to x 5 is only 0.008%. 

The x-dependence of the individual curves is rather flat, in particular close to x — 1 - 
except for z^ AM . The latter varies from +13.10 for x = to —1.52 for x = 1 with a 
zero for x ~ 0.81. It is interesting to note that the coefficient of C\Tn m is positive for 
< x < 1 whereas the contributions originating from the diagrams involving two closed 
fermion loops are negative and furthermore numerically smaller in a large part of the x 
interval. 

The boundary terms for x = and x = 1 can be extracted from the known single-scale 
three-loop results [9-11,13] and are given by 



,™(0) 


_ 7 FFL _ 


0.842 , 


FFM , 


;i) 


7 FFH _ 


2.494, 


^ AM (o) 


FAL ^ 


13.099, 


7 FAM, 


;i) 


FAH ^ 
~ 


-1.522, 


O) 


- 2z FLL ft 


i -3.916, 


FLM , 


;i) 


_ 7 FHL _ 


-0.067, 




__ 7 FHL ^ 


-0.067, 


FHM, 
m > 


;i) 


_ 9 FHH r 


a -0.384 


z™ M (0) 


FLL ^ 


-1.958, 


FMM , 
m ' 


;i) 


_ 7 FHH _ 


-0.192. 



(17) 



These limits constitute an important cross check of our calculation. We find perfect 
agreement with the known results. 

We refrain from listing explicit results for the x- 
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Figure 3: z^ FM , z^ AM , z^ LM , z^ HM and z™ IM as a function of a;. The exact curves are 
represented by the thick and the linear approximations for small x by the thin lines. 



dependent coefficients which can 



be 



obtained 



at 



the 



URL 

http : //www-ttp . physik . uni-karlsruhe . de/Progdata/ttp07/ttp07-21/ Instead we 
provide approximation formulae which have an accuracy of better than 1% for < x < 1. 
They are inspired by the expansion for small values of x and the behaviour for a? = lJU 
We obtain 



~FFM 
m 



~FAM 

Z<rn 



~FLM 



~FHM 

Z m 



0.842 + 4.333a; - 1.365s 2 + 3.136a; 2 log a; - 1.316a; 3 

13 3 

+L U (-0.972 + 1.820a; - 0.811a; 2 + 0.279a; 3 ) L 2 L 3 , 

32 32 

13.099 - 12.945a; + 9.041a; log a; - 1.676a; 2 

373 11 

+L, (6.724 - 4.407a; + 1.807a; 2 - 0.549a; 3 ) + — L 2 + —L\ , 



-3.916 + 2.948a; - 3.304a; log a; + 0.901a; 2 

13 1 

+L U (-1.921 + 1.604a; - 0.660a; 2 + 0.201a; 3 ) L 2 L 

v ' 36 ^ 18 • 

-0.067 - 0.612a; 2 + 0.443a; 3 - 0.148a; 4 

13 1 

+L U (-0.776 + 1.597s - 0.631a; 2 + 0.179a; 3 ) - — L 2 - —L 
v ' 36 ^ 18 ' 



3 



3 

H i 



3 We want to mention that the formulae provided on the web page are valid for < x < 5. 
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-FMM 



-1.958 + 0.501a; - 3.403a; logx + 1.264a; 2 - 0.103a; 3 log a; 

13 
72" 



13 1 

+L fJ- (-0.960 + 1.600a; - 0.653a; 2 + 0.198x 3 ) - — L\ - — L\ . (18) 



Let us in the following consider the result in the case of QCD, i.e., we set Cf = 4/3, Ca 
3, T F = 1/2, = 1, n m = 1 and define the quantities 



(2) ,M 
m 

(3) ,M 



C pTpz. 



FM 



(19) 



FFM n /i T FAM , ^ ^2 FML , ^ t^2 FMH , ^ t^2 FMM 
{ ^F 1 F z m ~v ^ F^ A 1 F Z m -\-L^F 1 F Z m -\- { ^F 1 F Z m ~T~ ( - y F 1 F Z m 



Evaluating the coefficients for the ^-expansion in numerical form the results become very 
compact and are given by (/i 2 = M 2 ) 



,(2),M 



,(3),M 



+ x 4 
+ x 8 
+ x c 
+ X 1 
+ x 2 
+ x 3 
+ x 4 



+ x° 
+ x 6 
+ x 7 
+ x 8 



1.0414) + x 1 (-1.6449) + x 2 (1.0000) + x 3 (-1.6449) 
1.2474 - 0.7222L* + 0.3333L 2 ) + a; 6 (-0.0844 + 0.0889^) 
-0.0118 + 0.0214^) , 
26.2712 - 1.3054^) 

-21.0921 + 16.9977L X + 1.0385^ - 1.0966^) 

12.7021 + 10.6870L X - 0.2222n { ) 

-13.0084 + 16.5103^. - 0.2408n, - 1.0966L x .n,) 

4.1035 + 0.1938L X - 0.2593L 3 + 0.7919L 2 + m (0.3613 - 0.2822L X . 
+ 0.074LL 3 - 0.2407L 2 )] 
-1.4908 + 2.8512L X + 0.4935^) 

0.1654 - 0.5756^ + 0.8224L 2 - 0.187372/ + 0.0267L x m) 
-0.1751 + 0.5452L X + 0.0653raj) 

-0.0705 + 0.0492L,, + 0.3125L 2 + n t (-0.0377 + 0.0025L X )] , (20) 



where the contributions proportional to ni are listed separately and L x = log x. In Fig. H] 
the expansions of the quantity Zm' M up to x n (n = 1, 3, 5, 8) are shown together with the 
exact expressions where n\ = 3 has been chosen corresponding to the case Mf = M c and 
M&. One observes a rapid convergence when including successively higher orders. 



M„ 



In Tab. [T] numerical results for the individual coefficient functions, for Zm M and Zm are 
shown in the region around x = 0.3 where again n\ = 3 has been adopted. 

Quite often it is convenient to consider the masses which only appear in closed loops in the 
MS scheme. Thus, transforming Mf to the MS scheme leads to the following modifications 
of Eq. (USD: 



Z rn' i X ) 



Z m' ( X f) J 

z^ M (x f )+C F T F n m AzUx f ) 



(21) 
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Figure 4: Zm as a function of x where n\ = 3 and jj 2 = M 2 has been chosen. The exact 
result is is shown together with the the expansions up to x n (n = 1,3,5,8) (dotted line 
to long-dashed line). 



X 


FFM 


FAM 


FLM 


FHM 


FMM 

m 


J3)M 


z {3) 


0.26 


1.57 


6.45 


-1.93 


-0.10 


-0.55 


12.15 


-111.66 


0.27 


1.59 


6.29 


-1.89 


-0.10 


-0.52 


11.89 


-111.92 


0.28 


1.61 


6.12 


-1.84 


-0.11 


-0.50 


11.62 


-112.19 


0.29 


1.63 


5.96 


-1.80 


-0.11 


-0.48 


11.37 


-112.44 


0.30 


1.64 


5.80 


-1.76 


-0.11 


-0.46 


11.11 


-112.70 


0.31 


1.66 


5.64 


-1.72 


-0.11 


-0.44 


10.86 


-112.95 


0.32 


1.68 


5.49 


-1.68 


-0.12 


-0.42 


10.61 


-113.20 


0.33 


1.70 


5.34 


-1.64 


-0.12 


-0.41 


10.37 


-113.44 


0.34 


1.72 


5.19 


-1.60 


-0.12 


-0.39 


10.13 


-113.68 



Table 1: Numerical results for x-dependent coefficients contributing to the three- loop 
quantity Zm . For the renormalization scale fj 2 = M 2 has been adopted. 



with 



12 



Az f m (x f ) = i (h^ + 2^ x{24x 3 log 2 (x) + I2x log(x) + 24x 



+ 7T 2 [x 2 (Ax - 9) - 3] - 6 [(4a; + 3)x 2 + l] [log(x) log(x + 1) + Li 2 (-x)] 
- 6{x - 1) (Ax 2 + x + 1) [log(l - x) log(z) + Li 2 (x)] } , (22) 

where l fJif = log(/z 2 /m 2 ). In Eq. fl22|) we introduced the scale /i/ for the renormalization 
point of the quark mass m/ which is different from \x implicitly present in Eq. Q. The 
latter contains the information about the running of a s whereas the former incorporates 
the anomalous mass dimension of rrif. 

Let us for completeness also present the inverse relation l/z m which is conveniently ex- 
pressed in terms of the MS quark mass m q . Using analogous conventions to ([19]) we 
find 



(l/z m )^ M (x q ) 
(V^) (3) ' M K) 



(M q /m q (^ M (x q 
(M q /m q (^ M (x q 



MM, 



X 



,(3),M 



{x q ) 



C 2 F T F n m Az q m {x 



(23) 



with 



X n 



Az* m (x q ) 



1 



ig2 {48(31, + 7)x q log 2 (x q ) + \44x\ \og(x q 

+ 8tt 2 [1x\ - 15x 3 - 3x q - 1) + 137 

3 



l, 



48 



312a; 



72x 2 q + 241, [-1, + 4 + 12vr 2 (-2x A q + 3x 3 q - 3x q + 2) + 13 



x n 



(7x q + 5)x 3 q + x q + -I, (2x1 + - , 
x [log(x,) log(x, + 1) + lA 2 (-x q ) 



Xn 



2) 



+ 4%{x q - 1) 



-I, (2x\ + x\ + x q + 2) + 7x 3 q + 2x\ + 2x q + 1 



log(l - x 9 ) log(x 9 ) - Li 2 (x g )] } . 



(24) 



It is understood that the renormalization scale dependent logarithms appearing in 
Zm' M {x q ), Zm' M (x q ) and Az%[x q ) are defined as I, = log(fi 2 /m 2 ). 

Again we can consider the masses which only appear in closed loops in the MS scheme, 
which leads to the following modifications of Eq. (T231 : 



(\jz n 
(l/z v 



)(2),M, 
v(3),M. 



X 



X n 



(l/z v 

(1/Zr, 



)(2),M i 
v(3),Af 



X fq)> 

(x fq ) + C F T F n m Az^(x fq ) 



(25) 
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with 



x fq = x fq {^») = 

Az f J{x fq ) = ^ + 2^ x fq { - 2Ax% \og 2 (x fq ) - 12x fq log{x fq ) - 2Ax fq 

+ 7T 2 [(9-4x fq )x 2 fq + 3] 

+ 6 [(4a;/, + 3)x 2 fq + l] [log(x fq ) \og(x fq + 1) + U 2 (-x fq )\ 

+ G(x fq - 1) [Ax% + x fq + 1] [log(l - x fq ) \og(x fq ) + U 2 (x fq )} } .(26) 

At the end of this Section we want to compare our result with the one of Ref. [14] where 
Zm has been computed in the linear approximation. Our result for the linear term 
reads 



(M f /m q (m q ))^ M (x q )\ linear = x q 



19.996 - 16.998 logx 9 + n t (-1.039 + 1.097 log x 9 

(27) 



We find agreement for three of the terms but the coefficient with the numerical value 
19.996 takes the value 21.277 in [14] @ This difference can be explained by the approxi- 
mations performed in Ref. [14] in order to extract the linear term of the mass relation^ 



4 Wave function renormalization constant 

In contrast to Z® s the wave function renormalization constant contains next to ultraviolet 
also infrared divergences. Thus it is not possible to construct a finite quantity by consid- 
ering the ratio between the on-shell and MS renormalization constant. For this reason we 
discuss in what follows the coefficients of the e-expansion separately. 



4.1 Two-loop result 



The two-loop corrections to Z 2 have been computed in Ref. [3]. We confirmed this 
result and obtain 



2 £\24 3 6 i 4 " Ue 36 3 



443 5tt 2 tt 2 7 2 5tt 2 , tt 2 , 

H 1 x H — x x 1 H x 

432 72 4 6 12 6 



4 Note that ni as introduced in Ref. [14] corresponds to our combination ni + n m . 
5 



We thank A. Hoang for communications on this point. 
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11 5 

- + -£ + - x 3 + x 4 ) [log a; log(l + x) + Li 2 (-x)] . {21 



( - + J logx+ ( 3 + ' i loL '" r 
11 5 

+ ( ~ o + o x + 7 x3 ~ %i ) [l°g x l°g(l — a; ) + Li 2 (a:)] 
3 2 6 

1 5 
+ -x + - 
3 2 6 

The convergence properties are very similar to Zm' M and shall not be discussed here. 
However, we would like to present a handy approximation formula which is obtained by 
an interpolation where the logarithmic divergence for x — > is extracted before. It reads 

Z? )M = -f---logxV 1 ^+f 1 + ---logx')LM 

2 e \24 3 / 4 M \6e 36 3 / 

+ - log 2 x + - logx + 1.711 - 2.356x + 1.125x 2 - 0.344a; 3 . (29) 
3 9 

and works to better than 1% for x E [0, 1]. 



4.2 Three-loop result 



We again refrain from listing explicit results for the x-dependent coefficients and refer 
to the URL http : / /www-ttp . physik . uni-karlsruhe . de/Progdata/ ttp07/ttp07 -2i~7| 
where the expressions can be downloaded in Mathematica format. It is, however, useful 
to present results for the analogue quantity to Zm as defined in Eq. (1T91 . The cubic 
and quadratic poles can be presented analytically and read 



r(3),M 



1-g 

96 



23 89 



108 



19 

— L a H loex + ni I 1 L u logx 

96 16 V 108 36 18 & 



1 

108 



( 1 1 r 1 i 

L„ H loex 

V 32 32 M 16 & 



(30) 



Concerning the single pole and the finite part we again present handy approximation 
formulae which we obtain by an interpolation to our expression after subtracting the 
singular terms for x — >• 0. We cast the result in the form 



,(3),Af 



r(3),M 



a + ai L^ + a 2 L\ + a 3 L 3 + £ (b + b x L^ + 6 2 L 2 + 6 3 Lj) , 

c + c x L^ + c 2 L 2 + e (d + + tfe-Lj) , (31) 
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where the coefficients are given by 



a 

ai 

a 2 

a 3 
bo 

bi 

b 2 

h 
c 

Cl 

c 2 
do 
d l 
d 2 



0.38426 log(x) [log(x) + 0.53600] [log(x) + 13.51219] 
+ 25.383 - 22.326a; + 11.127a; log(rr) - 1.473a; 2 , 

log(x) [-1.46528 log(a;) - 5.54630] + 6.832 - 1.971a: + 0.982a; 2 - 0.324a; 3 , 
403 , , , 127 

288 hg{x) + 48 ' 
193 

~576 ' 
407 7T 2 

864 + 128 



703) 
96 
9 



9 



— — — _ "i. los(x) — — log 2 (a;) , 

96 128 16 &V ; 16 & V ; ' 
9 9 

64 32 BK !l 
^ 3 

_ 64' 

log(a;) [-1.34028 log(x) - 2.41667] - 1.352 + 2.367a; - 1.180a; 2 + 0.387a; 3 , 
23 257 , , . 

36 + 144 l ° g{x) ' 
41 

~64' 

35 7r 3 . . 3 2/ \ 
-288 - 384 -16 l ° g{x) ~ 16 bg {x) ' 

I + ll l0g(x) ' 



J3_ 

64 



(32) 



For the singular contributions, which we know to high precision, we provide five digits 
after the decimal point whereas for the results from the fit three digits are given. 

As in the case of z m we want to present the expanded results for x — > which are given 
by (/i 2 = M 2 ) 



„(2),M 



,(3),M 



1 1 



. 443 5tt 2 4 , . . 2 , 2 . . 
- o lo g(^) ) + 4^ + -72" + 9 log ^ + 3 log ^ 



24 3 

— x + 2x 2 
4 



57T 2 

12 



125 7T 



11 



x 3 + — + — - — log(a;) + log 2 (a;) x 



,.22 16 , . v , 
F(-^ + ^log(a;))a; 6 + 



72 6 
379 



7840 28 



6 



+ — log(x) ) x 8 , 



1 



0.0104- 0.0104^ 



(33) 
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1 

+ - 

e 2 

1 " 
+ - 

e 



- 0.213 + (0.0093 - 0.0556L x )7i, + 1.18751^ + (0.0625L X + 0.0313)£ 
1.3978 + (0.0151 + 0.07412^ + 0.0556l£) n x - 2.6389L X - 1.5069L* 



2\ 4 
1 X 



+ 2.4674x - 2x 2 + 4.1123x 3 + (-3.381 + 1.8333L X - L 2 X ) 
+ (0.2933 - 0.3556L x )x 6 + (0.0483 - 0.1071L x )x 

+ (-0.1472 - 0.1875L X . - 0.1875Z,*) £ 

+ 31.2973 + (-1.9715 + 0.1765L X . - 0.0741L* + 0.037L 3 ) n x 
+ 2.25342,^ + 5.62042^ + 0.273124 



+ 



14.695 + 19.328L X + (1.5578 - 1.6449LJnj 



x 



+ (24.2836 - 0.7778n, + 27.2951^ - 0.63892^) 



2\ 2 



+ 
+ 

+ 



23.8364 + 28.4056L X + (-2.74162^ - 0.6021)nj 



x 



- 12.2965 + 7.3074L X - 2.4541L* + 1.5787L 3 . + (0.527 - 0.8466L, 



- 0.72222^ + 0.2222L 3 ) n x 



x A + (1.0836 + 1.7272n; + 8.6633L Jx 5 



- 2.8058 - 0.13422^ + 3.1252^ + 0.138924 + (-0.7078 + 0.0772^ 



x 



+ (0.4472 + 0.2937n; + 1.6073LJX 7 + 



- 0.9511 + 0.87872^ + 1.4563L 



x 



+ 0.099524 + (-0.1831 + 0.00542^.)^ 
+ (0.4605 + 0.8834L X + 0.562524 + 0.37524) £ . 



(34) 



In Fig. [5] we compare the exact results for the individual colour structures (for £ = 0) 
with the approximations including terms up to order x 2 . Only for x > 0.2 a difference is 
visible. 

For completeness we show the analogue curves for the coefficient of the 1/e pole in Fig. [61 
It is interesting to mention that only Z FFM has a non-trivial x-dependence beyond the 
logarithmic divergences for x — * 0. 



5 Applications and conclusions 

As an application of our result we want to discuss the charm quark effects in the relations 
between the pole, the MS and the IS" quark mass. For illustration we use m^rrib) = 
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0.2 0.4 0.6 0.8 m 1 



Figure 5: Finite parts of Z[ ML , Z[ MH , Z[ FM , Z[ AM , Z[ MM for Feynman gauge (f = 0) 
and L M = as a function of x. The exact curves are represented by the thick and the 
small-a; approximations by the thin lines. 



4.2 GeV, m c (m c ) = 1.3 GeV, \i = m b anc@ 
on-shell and the MS mass leads to 



( 4 ) 



m&) = 0.2247. The relation between the 



M h 



4.2 + 0.401 + 0.199 + 0.0094 



0.145 + 0.0182 



GeV , (35) 



where the tree-level, one-, two- and three- loop results are shown separately. The contri- 
butions from the charm quark mass which vanish for m c —>■ are marked by a subscript 
m c . One observes that the higher order contributions are significant. In particular, the 
two-loop charm quark effects amount to 9 MeV and the three-loop ones to 18 MeV. A 
similar bad convergence is observed in the relation between the IS" mass [36] M b ls and the 
pole mass M b . For M b = 4.8 GeV, m c (m c ) = 1.3 GeV, fj, = M b and a { ^{M b ) = 0.2150 it 
is given by 



M b is 



4.8 - 0.049 - 0.073 + 0.0041 



0.098 + 0.0112 



GeV. (36) 



6 As a starting point we use ots(Mz) = 0.118 and perform the running and decoupling with the 
program RunDec [35]. As in Ref. [36] we consider the mass relations with four active flavours. 
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Figure 6: Same as in Fig. [5] but for the \ je pole. 



However, the relation between the 15* and the MS quark mass is much better behaved as 
can be seen in the following example where we have chosen M b is = 4.69 GeV, m c (m c ) = 
1.3 GeV, /x = M b is and ai 4) (M fe ls ) = 0.2167 



4.69 - 0.382 



0.098 + 0.0047 



- 0.030 + 0.0051 



GeV. (37) 



The two-loop charm effects amount to only 4.7 MeV and three-loop ones to 5.1 MeV. We 
want to mention that in case only the linear approximation [14] of the charm quark mass 
effects is used the corresponding three-loop results in Eqs. (|35|) and (!37|) read 0.0167 and 
0.0037, respectively. 

We are now in the position to compare with Eq. (168) of Ref. [14] which provides the 



is 



relation between the 15* and MS bottom quark mass allowing for a variation of M b 
a s and the renormalization scale /i. Updating the coefficient of the m c term one obtains 
the formula 

m b = 4.169 GeV - 0.009 (m c (m c ) - 1.4 GeV) (38) 
which has an accuracy of better than 0.01% for 1.1 GeV < m c < 1.7 GeV. 
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A further application of our result would be the incorporation of our corrections in the 
analysis of the bottom quark mass determination from the T(IS') system. In the analysis 
performed in Ref. [37] the charm quark mass has not been considered and an uncertainty 
of ±10 MeV has been assigned which could be reduced to a large extend. 

To conclude, in this paper we have computed the three-loop QCD corrections to the 
on-shell renormalization constants for a heavy quark mass and the corresponding wave 
function where a second massive quark appears in closed loops. The two-scale three- 
loop diagrams are analytically reduced to 27 master integrals. The e-expansion of the 
latter is computed in analytical form, except for six coefficients for which one- and two- 
dimensional integral representations are available. We derived a compact expansion of 
the renormalization constants in the limit where the second quark mass is small. A rapid 
convergence is observed providing a good approximation to the exact result even close to 
the equal-mass case. 
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A Master integrals 



In this Appendix we collect the master integrals appearing in our calculation in graphical 
form. 




6.1 



s \ 





6.3 



Figure 7: Three-loop master integrals with six lines. Solid and dashed lines denote massive 
lines with masses m q and m/, respectively. Wavy lines are massless scalar propagators. 
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Figure 8: Three-loop master integrals with five lines. 
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